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We present an extensive experimental and theoretical study of surface acoustic wave-driven fer- 
romagnetic resonance. In a first modeling approach based on the Landau-Lifshitz-Gilbert equation, 
we derive expressions for the magnetization dynamics upon magnetoelastic driving that are used 
to calculate the absorbed microwave power upon magnetic resonance as well as the spin current 
density generated by the precessing magnetization in the vicinity of a ferromagnet/normal metal in- 
terface. In a second modeling approach, we deal with the backaction of the magnetization dynamics 
on the elastic wave by solving the elastic wave equation and the Landau-Lifshitz-Gilbert equation 
selfconsistently, obtaining analytical solutions for the acoustic wave phase shift and attenuation. 
We compare both modeling approaches with the complex forward transmission of a LiNbOa/Ni 
surface acoustic wave hybrid device recorded experimentally as a function of the external magnetic 
field orientation and magnitude, rotating the field within three different planes and employing three 
different surface acoustic wave frequencies. We find quantitative agreement of the experimentally 
observed power absorption and surface acoustic wave phase shift with our modeling predictions 
using one set of parameters for all field configurations and frequencies. 

PACS numbers: 76.50.+g, 75.30.Gw, 75.78.-n,75.80.+q 
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I. INTRODUCTION 

Spin mechanics exploits the coupling of spin degrees of 
freedom with elastic properties in order to control mag- 
netic properties and has been established in various ma- 
terial systems in the static and dynamic regime In 
magnetostrictive materials, a change in the magnetiza- 
tion orientation results in elastic strain and vice versa, 
enabling strain-controlled magnetization switchingi 6 -^ 
Furthermore, an acoustic wave propagating through a 
ferromagnet generates elastic strains allowing to drive 
magnetization dynamics via the magnetoelastic inter- 
action if the acoustic wave frequency is in resonance 
with the spin system at a given effective magnetic field. 
This magnon-phonon interaction has been investigated 
theoretically-^— and experimentall y 18 ' 19 in particular 
for bulk acoustic waves. 

Radio frequency (rf) elastic strains in ferromagnetic 
thin films can be generated by a surface acoustic wave 
(SAW) if the ferromagnetic film is deposited directly on 
the surface of a SAW-carrying crystal, resulting in a rigid 
elastic coupling. The physics of such heterostructures has 
been investigated experimentally^Sr— Recently, it has 
been demonstrated that the observed changes of the com- 
plex magnetotransmission through a hybrid LiNbOa/Ni 
SAW delay line can be identified with absorption and dis- 
persion signals of a SAW-driven ferromagnetic resonance 
(FMR) We here use the term FMR synonymously with 
spin wave resonance (SWR), which is the more appropri- 
ate term if the ferromagnet is large compared with the 
acoustic wavelength and a spin wave mode is excited. 

Since SAW-based devices are well established 2 ^ and 



can be fabricated by lithographic techniques they are at- 
tractive candidates for acoustically driven FMR in combi- 
nation with ferromagnetic thin films. In acoustic FMR, 
the magnetization dynamics is excited by a purely in- 
ternal magnetoelastic effective driving field, due to the 
rf strains generated by the acoustic wave in the ferro- 
magnet. The speed of sound is 5 orders of magnitude 
smaller than the speed of light, leading to a wave length 
of the SAW of the order of microns at GHz frequen- 
cies. This potentially allows for FMR experiments with 
micron-scale spatial resolution. Furthermore, spurious 
electromagnetic cross talk can be separated from mag- 
netoelastic effects associated with the acoustic wave by 
time-domain techniques . 24 ' 26 

The latter feature is of topical interest in the context of 
spin pumping, i.e., the generation of a spin current by a 
precessing magnetization at a normal metal/ferromagnet 
interface that can be detected by the inverse spin Hall ef- 
fect in the normal metal, e.g., a Pt layer j 2 ^— Since any 
ac electric fields which may be present at the position of 
the ferromagnet in a conventional FMR cavity can ham- 
per the interpretation of the spin-pumping signal^ 9 - - — 
it is important to separate electric and magnetic fields, 
which is naturally achieved by magnetoelastic driving as 
demonstrated recentl y 32 ' 33 To date, however, a thorough 
theoretical modeling of SAW-based FMR in combination 
with an extensive experimental study is still missing. 

In this work, we present a theoretical framework for 
SAW-FMR based experiments starting from the Landau- 
Lifshitz-Gilbert (LLG) equation and the elastic wave 
equation, providing expressions that can be used to ana- 
lyze SAW FMR and acoustic spin pumping experiments 
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and are applicable to various types of SAWs (Sec. HI1) . 
In a first approach, referred to as the "effective field ap- 
proach" , we show that a magnetoelastic driving field can 
be formally treated equivalently to a conventional, ex- 
ternal "tickle" field (Sec. Ill Ap . An expression for the 
magnetization dynamics upon magnetoelastic driving is 
found, with which it is straightforward to account for 
the particular type of SAW employed by including the 
relevant strain tensor components in the magnetoelas- 
tic contribution to the free enthalpy density. Based on 
this expression, the power absorbed in FMR is calculated 
f Sec. Ill A~Tj) , as well as the spin current generated by the 
magnetization precession at a ferromagnet / normal metal 
interface (Sec. Ill A"2]) . In a second approach, referred to 
as the "backaction approach" , we solve the LLG equation 
and the elastic wave equation simultaneously to obtain a 
intuitive physical picture of the backaction of the FMR 
on the elastic wave (Sec. IIIBI) . With the obtained ana- 
lytical solutions, both phase shift and attenuation of the 
elastic wave can be calculated. In addition, we present 
an extensive experimental study of SAW-FMR where the 
external magnetic field is rotated within three different 
planes employing three different SAW frequencies. We 
compare the experimental results with simulations based 
on the "effective field approach" and the "backaction 
approach" to discuss the advantage of each particular 
method (Sec. IIII[) . We demonstrate that the experimen- 
tal data can be modeled using a single set of simulation 
parameters for all measurement configurations and both 
simulation methods. Finally, we summarize our results 
and give an outlook on to further experimental and the- 
oretical investigations (Sec. IIV[) . 



II. THEORETICAL CONSIDERATIONS 

In this section, we provide the theoretical framework 
for the SAW-FMR experiments presented in Section 
Mil In the experiment discussed^ the SAW propagates 
through a delay line with a ferromagnetic film deposited 
between two interdigital transducers (IDTs) as schemat- 
ically depicted in Fig. Q] and explained in more detail in 
Sec. Mil We thus consider a two-layer system, consist- 
ing of a piezoelectric, single crystalline substrate and a 
poly crystalline, metallic ferromagnet. Rather than treat- 
ing this problem numerically, we will make the following 
physical assumptions and simplifications to obtain ana- 
lytical expressions for the acoustic wave attenuation and 
phase shift upon FMR as discussed in the following. The 
SAW penetrates the piezoelectric by a length 8^ of the 
order of the acoustic wavelength, c.f. Fig. [TJ which is 
typically a few microns at GHz SAW frequencies^ The 
thickness d = 50 nm of the ferromagnetic thin film on 
the other hand is much smaller than the acoustic wave- 
length. Therefore, we assume the elastic strains within 
the ferromagnetic film to be homogenous within the yz- 
pane (cf. Fig.[T|), thus treating the acoustic wave travers- 
ing the ferromagnet as a bulk acoustic wave, with strain 



components identical to the ones at the surface (z = 0) 
of the piezoelectric substrate. We focus on the acous- 
tic wave propagation and magnetization dynamics ex- 
clusively within the ferromagnet, rather than solving the 
acoustic wave equation and LLG equation for a two-layer 
system. We thus consider a ferromagnetic slab as de- 




FIG. 1: (a) Schematic illustration of the SAW- 
delayline/ferromagnet hybrid device used in the experiment. 
The aluminum interdigital transducers (IDTs) launch and 
detect a SAW, which traverses the delay line, periodically 
straining the ferromagnetic thin film, (b) The gray shading 
illustrates the SAW penetrating the piezoelectric by a length 
S, which is of the order of the acoustic wavelength. The 
thickness d of the ferromagnet is much smaller than the 
wavelength, justifying the assumption that the elastic strain 
within the ferromagnet is homogenous across its thickness. 

picted in Fig. [5] and employ a right-handed coordinate 
system (x,y,z), where x is the propagation direction of 
the SAW, y a transverse direction, and z the film normal. 

In the absence of further perturbations, the direction of 
the magnetization m = M /M is determined in a macro 
spin approach 3 — by the static free-enthalpy density (nor- 
malized to the saturation magnetization M) of the nickel 
film, given by^ 

G = —HoH ■ m + B^ml + B u (m ■ it) 2 — [ioH cx ■ m, (1) 

where H denotes an externally applied magnetic field, /iq 
is the vacuum permeability, = [1qM/2 represents the 
shape anisotropy of the thin film, B u is an anisotropy pa- 
rameter defining a uniaxial in-plane anisotropy along the 
unit vector u, and m x , m y , and m 2 are the components of 
the unit vector m; /j,qH cx = D s Am is the exchange field 
with the exchange stiffness D s and A = d x 2 + d y 2 + d z 2 is 
the Laplacian operator with respect to the spatial vari- 
ables x, y, z. 

In harmonic approximation, the elastic energy density 
W is given by-i^^ 

W = -CijkiSijEki, (2) 
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FIG. 2: Relation between the coordinate systems employed. 
The (x,y,z) frame of reference is spanned by the propaga- 
tion direction of the SAW, the transverse in-plane direction, 
and the normal of the ferromagnetic film. To solve the LLG 
equation, the (1,2,3) coordinate system is employed, with its 
3-direction corresponding to the equilibrium magnetization 
orientation. The inset shows the precession cone of the mag- 
netization, with the transverse magnetization components mi 
and m2. 



where we make use of the Einstein summation conven- 
tion. Cijki denotes the tensor of the elastic constants and 

£ ij = (ff~ + 7S"")/ 2 witn e { x 'y' z } are tne strain 
tensor components where Ui are the components of the 
mechanical displacement field. The dynamic, magnetoe- 
lastic contribution to the free-enthalpy density reads 3 ^ 

G d = bi[e xx (x 1 t)m 2 x +e y y{x,t)m'l J +e zz (x,t)m 2 z } 



2b 2 [e xy (x, t)m x m y + s xz (x, t)m x m z 
e yz {x,t)m y m z ], 



(3) 



with the magnetoelastic coupling constants b\ and b 2 and 
the SAW-induced dynamic strain tensor components e. 
Equation ([3]) holds for cubic symmetry of the ferromag- 
netic layer; for a polycrystalline film it is further simpli- 
fied by b\ — b 2 ?^ Depending on the particular type of 
SAW (e.g., Rayleigh wave^ or Love wave^l), some of the 
strain tensor components derived from the mechanical 
displacement field components 3 ^, are zero as discussed in 
the following. The total free enthalpy density normalized 
to the saturation magnetization of the film thus reads as 



gtot = Q + Qd + W fM. 



(4) 



Starting from this expression, we will in the following 
sections solve the LLG equation to obtain expressions 
for the magnetization dynamics and its backaction on 
the elastic wave in ferromagnetic resonance. 



A. Landau-Lifshitz-Gilbert approach 

The equation of motion for the magnetization direction 
m under the influence of an effective magnetic field i? e g 
is the LLG equatio n 38 i 39 



where 7 and a are the gyromagnetic ratio and a phe- 
nomenological damping parameter, respectively. In equi- 
librium, the magnetization orientation is parametrized 
by the polar angles 0o and 4>o for which the static free 
enthalpy density G [Eq. (TT])] is minimal. We introduce 
a new carthesian frame of reference (1,2,3) in which m 
points along the 3-axis and the 2-axis is in the film plane 
as shown in Fig. [5J The transformation matrix relating 
the two coordinate systems depicted in Fig. [2] is given in 
Appendix [A] 

Allowing for small deviations from the magnetization 
equilibrium direction mo, we write for the magnetization 
direction 




0{m\,m%) 



(6) 



where m,\jm% <C 1 in the (1,2,3) coordinate system. The 
effective magnetic field is given by 



V m G' 



(7) 



where V m = (<9 m i, d m2 , d m3 ) is the vector differential 
operator with respect to the components of m. By defi- 
nition of the (1,2,3) coordinate system, the static effective 
field is parallel to the 3-direction in the equilibrium. Fol- 
lowing the ansatz of Baselgia et al.^ we expand V m G 
at the equilibrium position of m, considering terms up 
to the first order in mi and m 2 . For the dynamic compo- 
nent of the effective field, we consider only terms, which 
are of zeroth order in m\ and m 2 . We thus find for the 
total effective field 



Gnmi + G 12 m 2 
fi H cH = - ( Gi 2 mi + G 22 m 2 
G 3 



Gf 
G* 
G d 3 



(8) 



and Gij — 



with the abbreviations Gi — d mi G\ m=mo 
dm i d Tnj G\ rn=mo and the required explicit expressions 
for these derivatives given in Appendix [A] Making a 
plane-wave ansatz for the transverse magnetization = 
m® exjp[i(kx — ut)] (i = 1,2), w being the angular fre- 
quency, k the wave number, and considering only the 
transverse magnetization components, we find for the 
LLG 



11— <JT3 — 



G12 



1 



12 + — 




Mo 



hi 
h 2 



(9) 



dtin — —7m x hqHqs + am x dtm, 



(5) 



Hohi = —Gf — — d mi G \m=m are the components of 
the effective driving field 

fM>hi= - 2bi sin O cos O cos 2 (f> + e vv sin 2 4> Q - e zz ] 

- 2b 2 [(£ X z cos (j) + £ yz sin0o)cos(20 o ) 

+ 2e xy sin O cos O sin (j) cos cj>o] (10) 
= + 26i sin O sin O cos cj> [e xx - e yy ] 

- 2b 2 [cos O (s yz cos 0o - £xz sin <f> ) 

+ e xy sm0 o cos(2(j) Q )]. (11) 
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In the derivation of Eq. ©, we have neglected terms 
quadratic in rrii and products of m, with the effective 
driving field; the latter are assumed to be of the same or- 
der as m 2 . Equation © can be solved for the transverse 
magnetization, resulting in 



(a) IP 




(12) 



12 

Gn — G3 — 




with 



D = ( Gu - G3 -^)(G 22 -G 3 -^) 

7 7 

2 



(13) 



In Eq. p2|) . x is the Polder susceptibility tensor de- 
scribing the magnetic response of the ferromagnet to 
small time-varying magnetic fields perpendicularly ori- 
ented to mo- x is a function of the derivatives of the 
static component of the free enthalpy density Eq. ^ 
while the driving field is determined by the dynamic con- 
tribution to the free enthalpy density Eq. In a con- 
ventional FMR experiment, the time-varying magnetic 
"tickle" field is provided, e.g., by standing electromag- 
netic waves in a microwave cavity and is oriented perpen- 
dicularly to the external magnetic field. As long as the 
magnetization is parallel to the static external magnetic 
field, which is the case if the external magnetic field is 
large compared to the anisotropy fields, the driving-field 
is always perpendicular to the magnetization, indepen- 
dently of the magnetization orientation. This is in stark 
contrast to SAW-driven FMR experiments, where the 
driving field components h\ and hi exhibit a pronounced 
m-dependence, characteristic of the type of the SAW, 
i.e., of which strain component dominates the dynamic 
free-energy contribution Eq. 24 A Rayleigh wave con- 



tains the strain components e xx , e xz , and e z 
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while 



a Love wave is a horizontally polarized shear wave with 
the dominant strain component e xy £^ For the Rayleigh 
wave, at the surface the components e xx and e xz are 
phase shifted by 90 deg. We will now discuss this in- 
dependence of the driving field exemplarily for the strain 
components e xx , e xy , and e xz separately, setting all other 
strain components equal to zero. To this end, we calcu- 
late the driving field n^h for each of the three geometries 
depicted in Fig. [3] that were used for our measurements. 
If m lies within the film plane (i.e. 6q = tt/2), referred to 
as IP configuration, we obtain the m-dependent driving- 
field components given in the appendix IA1 by Eqs. (|A7|I - 
(|A9|) . For m in the aiz-plane (i.e. 4>o — 0), referred to as 
OOP1 configuration, we find the expressions (|A1Q[) - (| A12|) 
and for the other out-of-plane configuration investigated 
(i.e. 0o = 7r /4), referred to as OOP2 configuration, the 
driving field components are given by Eqs. (IA13|) - (|A15I) . 

The magnitude of the driving fields for the different 
measurement configurations and strain components is 
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FIG. 3: Measurement geometries employed in the experiment, 
(a) In the in-plane configuration (IP) the external field is ap- 
plied within the film plane, (b) and (c) In the out-of-plane 
configurations OOP1 and OOP2, the field is aligned in planes 
perpendicular to the film plane as indicated in the figure. 



plotted in Fig. [4] as a function of the equilibrium orien- 
tation mo parametrized by (j>o and 80 . As Fig. @] demon- 
strates, the magnitude of the magnetoelastic driving field 
strongly depends on both the m-orientation and the 
dominant strain component associated with the partic- 
ular type of SAW. These characteristic fingerprints allow 
for a discrimination of the SAW-driven FMR from other 
driving mechanisms, such as electromagnetic crosstalk in 
free space, and for an identification of the relative ampli- 
tude of the strain components involved in the FMR exci- 
tation. As shown in Fig. @J for strain amplitudes of 10 -6 
and magnetoelastic constants of 25 T, typically found for 
a SAW delay line/Ni thin film hybrid, the magnitude of 
Hoh is 50 11T. Based on Eq. flT2"j) . we will in the following 
subsections quantify the power, which is absorbed by the 
ferromagnetic film upon resonance, and the spin current 
density, which is generated by the precessing magnetiza- 
tion at a ferromagnet /normal metal interface. 



1. Absorbed Radio-Frequency Power 

If the frequency of the SAW fulfills the ferromagnetic 
resonance condition, i.e., if the real part of the determi- 
nant [Eq. (|13p ] vanishes, the magnetoelastic interaction 
excites a resonantly enhanced magnetization precession, 
which in turn leads to a change in the SAW amplitude 
and phase. In other words, part of the SAW power is 
used to drive the magnetization dynamics. This is ex- 
pressed in terms of a complex change of the transmitted 
power AP, which we derive from Eq. (|12j) as 



AP 



2 



KiK ) x ( 



dV, 



(14) 



where Vq is the volume of the ferromagnetic film. 

The real and imaginary parts of AP are the disper- 
sion and absorption signals associated with the FMR, 
respectively^ Even though we have not considered an 
explicit backaction mechanism of the FMR on the SAW 
so far, Eq. (fT4"]) allows to quantitatively simulate the ab- 
sorbed radio frequency power Im(AP) in FMR, since by 
energy conservation, the power of the SAW has to be 
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FIG. 4: Polar plot of the driving field's magnitude \fj,oh\ for 
m in the film plane [panel (a)-(c)], out-of plane with (j>o — 
[panel (d)-(f)] and with <j>o = 7r/4 [panel (g)-(i)]. The distance 
from the origin indicates the magnitude of the driving field. 
The field components corresponding to the plots (a)-(i) were 
calculated by Eqs. (|A7|) - (|A15|) with fei = b 2 = 25 T and ey = 
10 -6 (see Sec. IIIip : the scale for the driving field applies to all 
panels. 



traversing the ferromagnet, various strain components 
with different amplitudes may superimpose, resulting in 
a more complex fioH dependence of the FMR signal. For 
a Rayleigh wave traversing an elastically isotropic mate- 
rial there are three non-vanishing strain components at 
the surface of the half space, namely e xx , e xz , and e zz , 
with a 7r/2 phase shift between e xx and and e xz £^ If the 
magnetization lies within the film plane, only the terms 
of the driving field containing e xx and e xz are relevant 
and superimpose with a phase factor of 7r/2; the result- 
ing SAW-FMR signal is shown in Fig.[5](b), revealing an 
asymmetry of the SAW-FMR pattern when compared to 
Fig. [5] (a). This asymmetry has been observed in SAW- 
FMR experiments ) 24 ' 33 and allows to experimentally de- 
termine the relative amplitudes of the strain components 
associated with the SAW. For the discussion of the out- 
of-plane orientations, we again focus on the longitudinal 
strain component and e xx . In the OOP1 configuration, 
the magnetization is oriented within the film plane along 
the x-axis for B u — 0, due to the large demagnetization 
field -Bd = 0.4 T (Sec. IIIII) . except for the case where 
the magnetic field points exactly along the z-axis. We 
thus expect the SAW-FMR signal to vanish for m par- 
allel to the x- and z-axes, c.f. Fig. H](d). F° r tue OOP2 
configuration however, the driving field only vanishes for 
(j,qH exactly parallel to the z-axis [cf. Fig. 2] (g)] and is 
nearly constant for other field orientations, since m is 
forced into the sample plane by the magnetic anisotropy, 
resulting in the SAW-FMR signal shown in Fig. [5] (c) . 



reduced by Im(AP) in resonance as will be shown in 
Sec. [TTTJ The phase change of the SAW upon FMR on 
the other hand, is assumed to be proportional to Re(AP) 
allowing for a qualitative discussion of the dispersion sig- 
nal based on Eq. (| 14[) . For a quantitative analysis of 
the phase shift of the SAW, the backaction of the mag- 
netization dynamics on the acoustic wave is considered 
explicitly in Sec. Ill Bl 

We will now qualitatively discuss the absorption and 
dispersion signals based on Eq. ([T4"]) for the IP and OOP2 
configuration. In Fig. [SJ AP is plotted as a function of 
the magnetic field orientation and magnitude with a spe- 
cific set of parameters given in the caption. Within the 
film plane, \ is independent of the field orientation for 
B u = 0. Assuming a purely longitudinal strain along x, 
for the IP configuration the signature of the SAW-driven 
FMR directly reflects the fourfold symmetry of the mag- 
netoelastic driving field [cf. Fig. @] (a)] with the SAW- 
FMR signal vanishing at ip = and if> = tt/2, c.f. Fig. [5] 
(a). For a pure in-plane shear strain e xy (not shown) the 
SAW-FMR signature is similar to the one in Fig. [5] (a) 
with the difference that the SAW-FMR signal vanishes 
for ip = 7r/4 and tjj = 3tt/4 [cf. Fig. g] (b)]. An out-of- 
plane shear strain e xz (not shown) results in a twofold 
symmetry of the SAW-FMR signal with the intensity 
vanishing at ip — tt/2. Depending on the type of SAW 
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FIG. 5: Normalized imaginary and real parts of the complex 
power AP calculated using Eq. (|14[1 for a pure strain along x 
{sxx 7^ 0) (a) and (c) and a superposition of e xx and e xz 
with e xz phase shifted by n/2 (b). The magnetic field is 
rotated within the film plane [(a) and (b)] and perpendicular 
to the film plane (c); the measurement configurations and 
angles ip axe defined in Fig. [3] The simulation parameters 
were B u = T, B d = 0.4 T, u = 2tt x 3 GHz, and a = 0.1. All 
other simulation parameters only affect the overall intensity 
of the SAW-FMR signal. 
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2. Spin-Current Generation 

At a ferromagnet/normal metal interface, the time- 
varying magnetization upon FMR generates an spin cur- 
rent density, which can be detected via the inverse spin 
Hall effect in the normal metal, e.g., in a Pt layer 
deposited on top of the ferromagnet i22r— It has been 
demonstrated that such a spin current can be driven 
acoustically by means of SAW-FMR. 33 Starting from 
Eq. (I12[) . we derive the spin current density expected at 
a normal metal/ferromagnet interface in the presence of 
magnetoelastic driving. The spin current density reads 



.27 



4tt 



Re(g^) 



dm 
~~dt 



with the spin mixing conductance g^£2. The dc com- 
ponent of the spin current is found by time-averaging 
Eq. CEP! 



^ C = ^ Re (s ri )Ini(Km 2 



mim|) , (16) 



which for a circular precession 771,2 = im\ simplifies to 



J 



dc 



f^Re(pT4-) sin r,— with the precession cone angle 
t. Equation (fToT) . together with Eq. (fT2"j) . allows to quan- 
tify the spin current density generated in SAW-driven 
FMR, where it is straight forward to account for the spe- 
cific type of SAW by superimposing the relevant strain 
tensor components when calculating the magnetoelastic 
driving field Eq. (fTTj) . 



Backaction of the ferromagnetic resonance on 
the acoustic wave 



In section III A[ we have discussed the SAW-driven 
FMR based on the LLG equation, which is the equa- 
tion of motion for the magnetization, and used a simple 
energy conservation argument to quantify the absorbed 
SAW power upon magnetic resonance. In this section, 
we will consider the propagation of elastic waves through 
the ferromagnet and will discuss the coupling of these 
waves to the magnetization dynamics. The equation of 
motion of an elastic displacement with components Uj is 
given by the elastic wave equation^ 3 - 



pd t 2Ui = d Xk a 



ik • 



(17) 



The elastic displacement components are given in the 
(x,y,z) coordinate system. In Eq. (TTTl) . p denotes the 
mass density, and a%k are the components of the stress 
tensor, which is derived from the elastic energy density 
byi2 ' 



dW 
de ik ' 



(18) 



In the absence of magnetoelastic interactions, the elastic 
wave equation thus readsi^ 



pd t 2Ui = Cikjr 



(19) 



In the presence of magnetoelastic interaction, however, 
G d M enters in Eq. (fT5| in addition to the elastic en- 
ergy density; note that Eqs. ([1]) and © have been nor- 
malized to M. Thus, the LLG and the elastic wave 
equation are coupled via the magneto-elastic interac- 
tion. In the following, we will solve this coupled sys- 
tem of equations by making a plane-wave ansatz with 
the wavevector along the x-direction. The magnetization 
and the elastic displacement are thus written as Mi t 2 = 
(15) Mj 2 exp [i(kx — ojt)] and u x . y . z = u xyz exp [i(kx — tot)] 



and the elastic mode equations linearized in M\ , M2 read 
as 



pu: u x = c\\k u x + 2ibiks'm8o cos0o x 
x [sin 0o -M 2 — cos 9 cos <t>oM{\ 



(20) 



P0J 2 Uy 



C44k 2 u y — 2i&2fcsin#o x 



x [2 sin 4>q cos 4>o cos 9qMi + cos 20oAf2] 
poj 2 u z — c^k 2 u z 

+ 2ikbi [sin </> cos 6qM% — cos 29 cos </>oMi] , 

where we have assumed cubic symmetry of the filmi^ 3 - For 
elastically isotropic media en = (AS 2 — ES)/(3S — E), 
and C44 = S, where E and 5* are Young's modulus and 
the shear modulus, respectively. 

Together with the linearized LLG from Eq. ©, we 
thus obtain a system of five coupled equations with the 
coefficients u x , u y , U z , M\, and M2. In order for these 
equations to be solved self-consistently, the determinant 
of the coefficients has to vanish yielding a polynomial 
equation in k. To obtain a simple physical picture of the 
backaction of the magnetization dynamics on the acoustic 
wave, we treat the acoustic modes Eqs. (|2U)) separately 
and neglect the interaction between these modes. 

(i) First, we consider a purely longitudinal acoustic 
wave, i.e. u y — u z = 0. By combining Eqs. (fT0|) . (|11|) . 
(HU), and HID, we find 



? / , Ft) 2 r O 



(21) 



+ X22W 2 - (X12 + X2i)wiw 2 }) k 2 ] u x = 



where we have introduced the abbreviations w\ — 
2 sin cos $o cos 2 4>q and W2 = 2 sin 9q cos tfio sin 0o ! v i — 
\fc\xfp is the sound velocity of the longitudinal wave 
in the absence of magnetoelastic coupling. Furthermoer, 
we have introduced a filling factor F < 1, which reduces 
the coupling of the modes, accounting for the fact that 
only a small fraction of the total volume traversed by the 
SAW is ferromagnetic, as illustrated in Fig. [TJ The cross- 
section A through which the acoustic power of the SAW 
is flowing is given by the penetration depth 5 of the SAW 
times the width w of the interdigital transducers. Since S 
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is of the order of the wavelength A, F = d/S will be of the 
order of d/X, with the thickness d of the ferromagnetic 
film. 

(ii) Conversely, we find for a purely transverse, in-plane 
shear wave (i.e., u x = 0, u z = 0) 



The power associated with the acoustic wave can be 
derived from its Pointing vector and reads as 3 - 

Pac = \a P v k uj 2 u x {x = 0) 2 |S 2 T n f , (26) 



K - - t 2 (i 



(22) 



+ X22W4 + (X12 + X2i)w 3 w 4 }) k 2 ] u v = 0, 

with the abbreviations W3 — 2 sin $o cos #0 sin 0o cos 4>o , 
W4 = sin 9q cos 20 q, and the sound velocity of the trans- 
verse wave vt = \J C44 / p in the absence of magnetoelastic 
coupling. 

(iii) For a pure out-of-plane shear wave (u x — 0, u y — 
0) we find 



F9h 2 

r 2 2/i r AU 2 r 2 
v t p p 



(23) 



+ X22w| - (Xl2 + X2i)w 5 w 6 }) k 2 ] u z = 0, 

with W5 = cos 0o cos 2#o and wq — cos #0 sin </>o . 

A Rayleigh type of SAW contains a longitudinal com- 
ponent and an out-of-plane shear component and prop- 
agates with the Rayleigh sound velocity i>r<2£ To obtain 
a simple physical picture, we will ignore the transverse 
component and assume the wave to be purely longitudi- 
nal with v\ = ur. This simplification is motivated by the 
observation that most characteristic features of the angle 
dependent SAW-FMR measurements can be simulated 
assuming a magnetoelastic driving field with a purely 
longitudinal strain.— 

Assuming that the changes of the mode due to mag- 
netoelastic interactions are small and neglecting the ex- 
change interaction \D S = in Eq. (|A4[) ]. we find for the 
wave number of the perturbed longitudinal wave Eq. (|21[) 



k 
Ak 



ko 



= F- 



h Ak 



{Xn w i + X22W2 - (x 



12 



(24) 

X21V1W2}), 



with the unperturbed wavenumber fco = w/vr. Note 
that at the largest SAW frequency employed in this work 
(2.24 GHz) the exchange term D s k 2 with D s = 2.1 x 
2Q-17 rp m (p_ e f_ [4]] ) results in an isotropic resonance 
field shift of only 0.4 mT justifying the disregard of D s . 

With Eq. ([24| we thus find for the longitudinal acoustic 
wave in the presence of magnetoelastic interaction 

Ux (x = x ) = u x {x = oy koXo exp (iAkxa), (25) 



where xq is the length of the ferromagnetic film. In 
Eq. (|25p. we have defined a normalized complex scat- 
tering parameter S??™ 11 , which describes the resonant 
attenuation and phase shift of the acoustic wave with 
S™ im = 1 off resonance, c.f. Eq. (JUJ). 



where P° c is the acoustic power of the wave out of res- 
onance. For small deviations of \S^ m \ 2 = e -2im(Afe)*o 
from unity, the change of the acoustic power in ferromag- 
netic resonance can be shown to be identical to the imag- 
inary part of Eq. (|14[) setting all strain components but 
e xx = iku x to zero. Accordingly, the real part of Eq. (fT4|) 
can be related to the phase arg(S , 2 1 >rm ) = Ke(Ak)xo. We 
thus find 



AP = -2Akx Q Pl c 



(27) 



demonstrating the consistency of the two models dis- 
cussed in Sees. Ill Al and iHBl We stress that this relation 
holds for small perturbations Afc of the wavenumber and 
for a purely longitudinal mode. Particularly, larger phase 
changes of the order of 7r/2 as observed in the experiment 
cannot be accounted for within this modeling approach 
because they would lead to Re(AP) > P° c . 



1. Magneto-phonon polaritons 

To illustrate the coupling of the elastic wave and the 
magnetization dynamics, we consider Eq. (|21[) with the 
magnetization oriented within the plane, i.e., 9q = tt/2; 
for simplicity, we assume £?<j = B u = 0, and a = 0. 
Equation (f2"Tj) simplifies to 



2 ; 2 



4b 2 sin 2 (f>o cos 2 4>oPqH PqM 



\ 



\ 



4 M op((poP) 2 -(^)' 



(28) 



which can readily be solved for uj. The magnetic field- 
dependent part of Eq. (|28|) is proportional to the square 
of the driving field for m within the plane and only 
e X x ^ 0, cf. Eq. (|A"7|) . Considering Fig. g] (a), we ex- 
pect the strongest interaction of the FMR mode with 
the elastic wave for </> = 7r/4 and no interaction for 
0o = and 0o = 7r /2- We assume a magnetic field of 
poH = 73.5 mT, corresponding to a resonance frequency 
of 2.24 GHz, i.e., the highest SAW frequency employed 
in this work, and use typical parameters for nickel with 
the references given below; ur = 3440 m/s, d = 50 nm, 
h = 23 T, p = 8900 kg/m 3 , 7 = 2.185 ^ B /h with Bohr's 
magneton /xb and the reduced Planck constant K. The 
wavelength of the acoustic wave at the expected FMR 
position is A rcs = 2ttvr/(poH^) = vr/(2.24 GHz) and 
we therefore assume a filling factor of F = d/X Tes - We 
solve Eq. (|28|) for w and plot in Fig. [5] the obtained re- 
sult against k for 0o = 0, 0o = 7r /4, and 0o = k/2. 
As long as the acoustic wave frequency is off-resonant 
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with the FMR frequency at the given field, we observe 
a linear dispersion uj = v^k of the acoustic wave and 
a fc-independent mode at uj = /iqH"/. Because there is 
no coupling of these modes for — and (fto = 7r/2, 
we observe a mode crossing. For cf>o — 7r/4, however, 
the coupling of the modes is strongest and we observe a 
mode hybridization, manifesting itself in an avoided level 
crossing, as shown in Fig. [B] these hybridized modes are 
referred to as magneto-phonon polaritons. 42 As it can be 
seen in Fig. [51 the mode we are referring to as FMR mode 
is fc-independent and exhibits the avoided level crossing 
with the acoustic wave mode at k ^ 0, which is why it is 
more correct to refer to this mode as SWR mode. 




k(10 4 cm" 1 ) 



FIG. 6: Dispersion relation of the acoustic wave and the FMR 
mode, calculated by solving Eq. (|28l) with the parameters 
given in the text. The dashed lines correspond to the FMR 
mode and the acoustic wave for <f>o = and <j>o = t/2, where 
the coupling term in Eq. (|28[) vanishes. The solid lines repre- 
sent the modes for 4>o = vr/4 where the coupling of the modes 
is strongest, manifesting itself in the avoided level crossing. 



Similar to the modeling of a Rayleigh wave by a longi- 
tudinal bulk mode discussed above, the attenuation and 
phase shift of other SAW modes in FMR can be modeled 
taking into account the transverse bulk modes. A SAW- 
FMR experiment with a Love wave, e.g., can be modeled 
by an in-plane shear wave Eq. (|22j) . 



III. EXPERIMENT AND DISCUSSION 

Having established the theoretical framework of SAW- 
FMR, we turn to the experiment. We investigate the 
hybrid SAW delay line device schematically depicted in 
Fig. [U fabricated from a y-cut z-propagation LiNbC>3 
substrate. The 70 nm-thick aluminum transducers with 
an interdigital spacing of 5 /zm and a metallization ra- 
tio of 50% were fabricated using optical lithography and 
e-beam evaporation. This geometry results in a funda- 
mental frequency of 172 MHz, corresponding to a sound 



velocity of = 3440 m/s^ and yields a bandpass at 
odd harmonic frequencies; for the SAW-FMR experi- 
ments, frequencies of 0.86 GHz, 1.55 GHz, and 2.24 GHz 
were employed. A polycrystalline Ni film with dimen- 
sions d = 50 nm, w = 400 /im, and xq = 570 /im was 
deposited between the transducers by e-beam evapora- 
tion. The complex forward transmission S21 of the de- 
lay line, defined as the voltage ratio of the electromag- 
netic wave detected and applied at IDT 2 and 1, respec- 
tively, was measured using vector network analysis with 
an input power of P — 0.1 mW. The SAW-transmission 
was isolated from spurious signals such as electromag- 
netic cross-talk and multiple transit signals by Fourier 
transformation and time gating i 24 ' 26 The delay line was 
mounted between the poles of a rotatable electromagnet 
and measurements were carried out in the three different 
measurement configurations shown in Fig. [3] All exper- 
iments were performed at room temperature. For each 
harmonic frequency, the power of the SAW was deter- 
mined by measuring the scattering parameter Su, de- 
fined as the voltage ratio of the reflected and applied 
electromagnetic wave at the transducer. The fraction of 
the applied electromagnetic power coupled into the SAW 
is given by P S aw = A|Su| 2 P /2, where A|Su| 2 is the 
amplitude of the dip occurring at the corresponding har- 
monic center frequency of the transducer and Pq is the 
input power applied to the transducer; the factor 1/2 ac- 
counts for the transducer bidirectionality. The change of 
phase and amplitude of the SAW upon ferromagnetic res- 
onance was determined by normalizing the scattering pa- 
rameter 5§f m = S 2 i(fioH)/S 2 i(fioH oS ), where the off- 
resonance field /Lto-Hoff = 150 mT for the in-plane mea- 
surements and /logoff = 1-2 T for the out-of-plane config- 
urations. The resonantly absorbed power was calculated 
from the normalized measured scattering parameter by 

Pabs = (1 - |5 2 T m | 2 )P S AW. 

We will now compare the experimental data obtained 
for the absorbed power P a bs with simulations based on 
the imaginary part of Eq. (fl"4")) , which has the advantage 
that different strain components can be easily incorpo- 
rated in the simulation. Furthermore, we will compare 
the measured normalized scattering parameter S^y™ 
with the one simulated with Eq. (p?5"]) . i.e., describing at- 
tenuation and phase shift of a purely longitudinal mode. 

We start with the in-plane configuration where /iqH 
was rotated within the film plane. Here, ip denotes the 
angle between the SAW propagation direction and fioH 
as defined in Fig. [3] For each employed frequency, the 
measured and simulated data are shown in a false color 
plot with the identical scale for simulation and experi- 
ment. The top row shows the absorbed power, simulated 
with Eq. HU), i.e., the "effective field approach". The 
best agreement between experiment and simulation for 
all frequencies was found for the parameters a = 0.1, 
B d = 400 mT, B u = 2.5 mT, and D s = 0; the exchange 
stiffness was neglected for reasons discussed above. The 
value for a is about a factor of two larger than liter- 
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FIG. 7: Experiment and simulation of the angle-dependent SAW-FMR in the IP configuration. The definition of the angle ip 
is given in Fig. [3] (a) . In the top row the power absorbed upon FMR is shown for all employed frequencies. For the simulation 
the imaginary part of Eq. (|14[1 was used with all parameters given in the text and in Tab. U In the second and third row, 
the measured magnitude and phase of the normalized scattering parameter SJi™ are shown together with the corresponding 
simulation, using Eq. (f25|) . 



ature values for polycrystalline nickel thin films4^ We 
conjecture that other, non-Gilbert type, line broadening 
mechanisms play a role. These mechanisms could include 
two-magnon processes^ non-uniform excitation due to 
the small acoustic wave length^, non-resolved standing 
spin-wave modes^£~— or excitation of other spin-wave 
modes and their damping. The clarification of this issue 
would require a systematic, frequency-dependent study 
of the SAW-FMR linewidth, particularly at higher SAW 
frequencies where frequency-independent contributions 
to the linewidth are less dominant and where standing 
spin-wave modes can be spectrally resolved. However, 
this would require the fabrication of higher frequency 
SAW devices by electron-beam lithography, which is be- 
yond the scope of this work. 

The magnitude of the strain components derived from 
the simulations depends on the frequency with the re- 
spective values given in Tab. HI for higher frequencies, 
a larger fraction of the electromagnetic power applied 
to the transducer is emitted into the free space and 
therefore less electromagnetic power is transfered into 
the SAW. As further simulation parameters, the fol- 
lowing literature values were used: M = 370 kA/m 
(Ref. S), h = b 2 = -3A s c 44 = 23 T (Ref. 35), with 
the isotropic magnetostriction constant A s = —38 x 10~ 6 
(Ref. 1 5 01) and the elastic shear module C44 = S = 74 GPa 
(Ref. p|, p = 8900kg/m 3 (Ref. H), and 7 = 2.185p, B /H 
(Ref. 1531 ). As shown by the graphs, the characteristic 
angle-dependence as well as the absolute value of the ab- 
sorbed power of the experimental data is reproduced by 
the simulation with one set of parameters. 



In the second and third row of Fig. the magnitude 
and phase of the measured normalized S^i™ parame- 
ter are shown together with the simulated S%° rm using 
Eq. (|25p . i.e., the "backaction approach". For this sim- 
ulation, the same parameters as above were used. The 
strain components do not enter in Eq. (|25p . but F is 
used as a free parameter, which is of the order of d/X 
and therefore depends on the frequency; the frequency- 
dependent parameters used for the simulations are also 
summarized in Tab. [I] In the simulation, F was chosen 
such that the simulated phase arg^Jf™ 1 ) quantitatively 
agrees with the experiment. In the SAW-FMR disper- 
sion, the characteristic angle-dependence is reproduced 
in the simulation. With this set of parameters, the simu- 
lated magnitude |5 , 2 1 i rm | is also in reasonably good agree- 
ment with the experiment, given the simplicity of the 
model^i Note that the nearly field-orientation indepen- 
dent change in the magnetotransmission observed in the 
absorption and dispersion data at low positve fields stems 
from hysteretic magnetization switching) 21 ' 2 and can be 
included in the modeling by searching for the local en- 
ergy minimum of the free enthalpy density Eq. ([1} with 
respect to the magnetization direction rather than the 
global minimum. However, we here focus on the non- 
hysteretic SAW-FMR signature and thus disregard the 
hysteretic switching in the simulation. 

Particularly at 2.24 GHz, the simulated maximum 
value for |S , 2 1 i rm | overestimates the acoustic wave attenu- 
ation observed in the experiment. Most likely, the agree- 
ment could be improved by modeling the Rayleigh wave 
more realistically by considering the attenuation and 
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phase shift of a coupled longitudinal and transverse wave. 
Since |S£f rm | substantially deviates from unity in the ex- 
periment, the approximation discussed in Sec. Ill Bl when 
showing that the change of acoustic power for the lon- 
gitudinal mode upon FMR is identical to the imaginary 
part of Eq. ([M)) is not justified. Nevertheless, the mod- 
eling of the absorbed power with Im(AP) can reproduce 
the experimentally observed lineshape of the absorption 
data better than the modeling with Eq. (|25|) . since all 
relevant strain components can easily be accounted for. 
As already mentioned in Sec. Ill BL the experimentally 
observed phase shift of more than 90 deg., however, can 
not be quantitatively modeled with Eq. (1T41) . We finally 
note that, within both approaches, for the magnetic field 
orientations close to the SAW propagation direction, i.e., 
for ip close to deg. in Fig [7J the agreement between 
experiment and simulation is less good than for other 
orientations. This could be a consequence of the simpli- 
fications involved in the modeling and might be improved 
by considering non-linear driving] 54 ' 55 

TABLE I: Frequency-dependent parameters used for the sim- 
ulations in Figs. [7J51 

0.86 GHz 1.55 GHz 2.24 GHz 



£x* (10~ b ) O Ll5 0.36 

e xg (10~ 6 ) 0.09i 0.0575i 0.018i 

e zz (10" 6 ) 0.18 0.115 0.036 

F (d/X) 0.26 0.35 0.38 

A (/xm) 20/5 20/9 20/13 



Figures [5] and [HI show the experimental and simulated 
data for the OOP1 and OOP2 field orientations, with the 
angle tp defined in Fig. [3] and the panels organized in the 
same way as in Fig. [7J For the simulations the identical 
parameters as for the IP configuration were used, result- 
ing again in a good agreement between experiment and 
simulation. For both configurations a slight asymmetry 
along the ip-axis is observed. It stems from a misalign- 
ment of the sample with respect to the magnetic field of 
less than one degree and is reproduced in the simulation. 

We finally note that the overall agreement between ex- 
periment and theory for the absorption data is slightly 
better in the "effective field approach" than in the "back- 
action approach", particularly the quantitative value of 
the absorbed power can be reproduced very well in this 
approach, while the "backaction approach" is better suit- 
able for describing the phase. 

IV. SUMMARY AND OUTLOOK 

In summary, we have provided a theoretical frame- 
work for SAW-based FMR experiments based on two 
different approaches. In the "effective field approach", 
we calculated the magnetization dynamics in the pres- 
ence of a magnetoelastic driving field. Based on the ob- 
tained analytical expression, the rf power absorption in 



the ferromagnetic film can be calculated under magnetic 
resonance conditions. The advantage of this method is 
that arbitrary types of SAWs or bulk waves which may 
drive the ferromagnetic resonance can be modeled in a 
rather simple fashion. The modeling is achieved by su- 
perimposing the corresponding strain components in the 
magnetoelastic terms of the free enthalpy density, where 
necessary with an additional phase factor. Further, we 
have derived an expression for the spin-current genera- 
tion upon FMR, which applies to SAW-based acoustic 
spin-pumping experiments in ferromagnet/normal metal 
hybrids, again with the capability of modeling arbitrary 
acoustic wave modes. 

In the "backaction approach" , we have taken into ac- 
count the backaction of the magnetization dynamics on 
the acoustic wave by solving the LLG equation and the 
bulk clastic wave equation, neglecting the interaction be- 
tween the longitudinal and transverse bulk waves. De- 
spite the approximations involved, the analytical results 
allow a quantitative analysis of the phase shift and at- 
tenuation of the SAW upon FMR and we showed that for 
small perturbations of the longitudinal elastic wave the 
two approaches are consistent. The particular advantage 
of this approach is that the SAW phase shift in FMR can 
be modeled quantitatively. 

Furthermore, we have performed systematic SAW- 
FMR experiments using a Ni/LiNb03 hybrid SAW de- 
layline and measuring the complex scattering parameter 
of the device S%i rm as a function of the magnetic field 
orientation and magnitude; the field was rotated in the 
in-plane configuration and in two different out-of-plane 
configurations for three different SAW frequencies each. 
We have shown that the absorbed power can be quanti- 
tatively described with the LLG approach with one pa- 
rameter set for all measurement configurations and fre- 
quencies. Moreover, we have shown that the phase shift 
of the SAW upon FMR can be quantitatively modeled 
by the backaction model, considering a longitudinal bulk 
wave traversing through the hybrid. 

This work thus lays theoretical foundations for SAW- 
based spin mechanics experiments such as SAW-FMR 24 
driven acoustic spin pumping 33 and is applicable to var- 
ious (surface) acoustic wave modes. Further theoretical 
work could be directed toward describing the SAW prop- 
agation through the coupled two-layer system consisting 
of the piezoelectric substrate and the ferromagnetic layer 
numerically, i.e., without the approximations discussed in 
Sec.[ITl to obtain an even more accurate description of the 
experiment and quantitative modeling of the SAW atten- 
uation and phase shift. Particularly, it would be desirable 
to improve the agreement between theory and experiment 
for field orientations close to the propagation direction 
of the SAW. A particularly appealing experiment would 
be to spectroscopically resolve the avoided level crossing 
of the SAW and FMR mode. To this end, a ferromagnet 
with a FMR linewidth of the order of 10 MHz with a sim- 
ilar magnetoelastic coupling as Ni or SAW devices with 
higher frequencies and thus larger filling factors would be 
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FIG. 8: Experiment and simulation of the angle-dependent SAW-FMR in the OOP1 configuration. The definition of the angle 
if) is given in Fig. 0(b). In the top row the power absorbed upon FMR is shown for all employed frequencies. For the simulation 
the imaginary part of Eq. (|14p was used with the same parameters as in the IP configuration. In the second and third row, 
the measured magnitude and phase of the normalized scattering parameter SJi™ are shown together with the corresponding 
simulation, using Eq. (|25[) . 
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FIG. 9: Experiment and simulation of the angle-dependent SAW-FMR in the OOP2 configuration. The definition of the angle 
ij) is given in Fig.[3](c). In the top row the power absorbed upon FMR is shown for all employed frequencies. For the simulation 
the imaginary part of Eq. ()14[) was used with the same parameters as in the IP configuration. In the second and third row, 
the measured magnitude and phase of the normalized scattering parameter SJf™ are shown together with the corresponding 
simulation, using Eq. (|25[) . 



required. Therefore, it could be beneficial to intention- 
ally use standing spin-wave resonances, e.g. by employing 
magnonic crystals, and couple these to a SAW resonator. 
Alternatively, the coupling of paramagnetic centers to 



elastic waves^ e.g. via a crystal field^ could be ex- 
ploited to observe a strong coupling of the spin ensemble 
to a SAW resonator, in analogy to the strong coupling of 
spin ensembles to photon cavities £2r— 
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Appendix A: Coordinate transformation and free 
enthalpy derivatives 

The transformation between the (x, y, z) coordinate 
system, defined by the propagation direction of the SAW 
x and the surface normal z, and the equilibrium system 
(1,2,3) is given by 



m x \ 




( mx 


m y 


-1 




m z 1 




\ m 3 



(Al) 



with 



U = 



COS 8q COS 0o 

cos 9q sin (j>o 
— sin 8q 



— sin 4>o 
cos 4>o 




sin 9q cos (po 
sin #o sin </>o 
cos O 



(A2) 



The driving fields derived from the derivatives of the 
static free-enthalpy density Eq. JT]) with respect to the 
magnetization components are given by 



G, 
G21 
G22 



^7713 G\m— rr i 

-Mo 



H 3 + 2B d cos 2 6 Q + 2B xl u\ 



G\2 — ^mi C^mo ^\m—mo 
)2 n\ on _ : „2 



d 2 mi G\ r 



(A3) 

D s k 2 , 
2B u u 2 ui, (A4) 
2B d sm 2 8 + 2B u u 2 , (A5) 



(A6) 



In the following, the derivatives of the dynamic free en- 
thalpy density used for the plots in Fig. [4] are stated ex- 
plicitly for the strain components e xx , e xy , and e xz sepa- 
rately, setting all other strain components equal to zero. 
If m lies within the film plane (i.e. 80 — ir/2), referred to 



as IP configuration, we obtain the m-dependent driving- 
field components 
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(A7) 
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and 
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(A9) 



2b 2 e xz cos <f> 

h 2 J V 
Conversely, if m is in the xz-plane (i.e. 4>q — 0), re- 
ferred to as OOP1 configuration, we find 
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(All) 



(A12) 



For the other out-of-plane configuration investigated 
(i.e. 4>o = 7r/4), referred to as OOP2 configuration, the 
driving fields read as 



and 
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